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We calculate the most massive object in the Universe, finding it to be a cluster of galaxies with 
total mass M200 = 3.8 x 10 15 Mq at z = 0.22, with the la marginalized regions being 3.3 x 10 15 Mq < 
M200 < 4.4 x 10 15 Mq and 0.12 < z < 0.36. We restrict ourselves to self-gravitating bound objects, 
and base our results on halo mass functions derived from N-body simulations. Since we consider 
the very highest mass objects, the number of candidates is expected to be small, and therefore each 
candidate can be extensively observed and characterized. If objects are found with excessively large 
masses, or insufficient objects are found near the maximum expected mass, this would be a strong 
indication of the failure of ACDM. The expected range of the highest masses is very sensitive to 
redshift, providing an additional evolutionary probe of ACDM. We find that the three most massive 
clusters in the recent SPT 178 deg 2 catalog match predictions, while XMMU J2235. 3-2557 is roughly 
3a inconsistent with ACDM. We discuss Abell 2163 and Abell 370 as candidates for the most massive 
cluster in the Universe, although uncertainties in their masses preclude definitive comparisons with 
theory. Our findings motivate further observations of the highest mass end of the mass function. 
Future surveys will explore larger volumes, and the most massive object in the Universe may be 
identified within the next decade. The mass distribution of the largest objects in the Universe is a 
potentially powerful test of ACDM, probing non-Gaussianity and the behavior of gravity on large 
scales. 



Introduction — Our Universe has a finite observable vol- 
ume, and therefore within our Universe there is a unique 
most massive object. This object will be a supcrcluster of 
galaxies. Theoretical studies of the growth of structure 
have now matured, and the mass of the most massive 
objects can be robustly predicted to the level of a few 
percent. Furthermore, we are in the midst of a revo- 
lution in our ability to conduct volume-limited samples 
of high-mass clusters, with Sunyaev-Zel'dovich (SZ) and 
X-ray surveys able to provide complete samples at mass 
> 5 x 10 14 Mq out to z > 1. The masses of the most 
massive clusters in the Universe are therefore a robust 
prediction of ACDM models, as well as a direct observ- 
able of our Universe. 

The cluster mass function is already being utilized as a 
probe of cosmology, and in particular, of the dark energy 
equation-of-state PHD]. What additional value is there 
in singling out the very tail end of the mass function, 
representing the most massive clusters in the Universe, 
for special treatment? First, we note that these systems 
are in many ways the easiest to find, as they are among 
the largest and brightest objects. They thus avoid many 
selection effects which might plague lower mass cuts. In 
addition, these systems constitute a very small sample 
(ideally, just one compelling candidate), and it is possible 
to devote significant observational resources to studying 
them. One might imagine coupled S-Z, X-ray, and weak 
lensing measurements, and thus the masses of these sys- 
tems will be among the best constrained of any systems. 
The mass-observable relation for clusters is an essential 
component in using the cluster mass function to mea- 
sure properties of the dark energy, and therefore there 
is a tremendous amount of ongoing work to characterize 
the masses of these objects d, I8l-fl3l|. Finally, because 
we are probing far down the exponential tail of the mass 
function, these objects offer an unusually powerful con- 
straint. If the most massive object is found to have too 
large a mass (or especially, as explained below, too small 



a mass), this single object will provide a strong indication 
of non-Gaussianity or modified gravity (TJ]. An excel- 
lent example of this is the high-redshift cluster XMMU 
J2235.3-2557 (hereafter XMM2235) [lj], which has been 
argued to be a few sigma inconsistent with ACDM [Tsl — 
1 1 81 ] - A similar approach based on strong lensing has 
been presented in [19(, which considers the distribution 
of the largest Einstein radii in the Universe as a probe of 
ACDM. Although much work has focused on using halo 
statistics as a probe of cosmology, here we focus on using 
the high-mass tails of precision mass functions to make 
explicit predictions for current and future observations. 

A critical question in one's attempt to determine the 
most massive object is to define precisely what is meant 
by "object". The largest structure in the Universe de- 
tected to date is the Sloan Great Wall [2(| , but the iden- 
tification of this wall as a unique object is sensitive to 
a (completely arbitrary) density threshold. For our pur- 
poses we define an object as a gravitationally self-bound, 
virialized mass aggregation. These objects have decou- 
pled from the Hubble flow, and represent large local mat- 
ter overdensitics. This definition has the convenience of 
robustly identifying objects (both in theory and observa- 
tion). 

Mass function — Recent years have shown tremendous 
progress in characterizing the mass function of dark mat- 
ter halos in cosmological N-body simulations. We have 
now established, to better than 5%, the expected number 
density of dark matter halos as a function of mass and 
redshift [2l]-[23| . In the simulations underlying these pre- 
cise mass function expressions, the halos at the high-mass 
end are resolved by millions of particles, lending particu- 
lar confidence and robustness to the mass function in this 
regime. The simulations are pure dark matter, and ne- 
glect the influence of baryons. At smaller scales baryons 
could play a major role in the density profile of the dark 
matter halos, and could potentially impact the mass func- 
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tion of the objects themselves. At the large scales being 
considered in this paper, the effects of baryons are ex- 
pected to be negligible. This is particularly true as our 
interest is in the mass function, and hence the number 
density of these halos, not their density profiles. 

An important issue is the process by which a dark mat- 
ter halo is identified and characterized in a dark mat- 
ter simulation [24j]. There are two dominant approaches: 
friends-of- friends (FOF) and spherical overdensity (SO). 
FOF defines a halo by contours of constant density, while 
SO defines halos by the overdensity (compared to the 
mean or critical density) within a spherical region. It 
has been argued that the mass associated with SO can 
be most closely tied to observations of clusters [23j . On 
the other hand, using an FOF with a linking length of 0.2 
corresponds closely to contours of density 200 times the 
background density, which from spherical collapse mod- 
els is a natural proxy for the virial mass. Because of the 
steep exponential in the mass function, our results are 
essentially independent of these differences (see Fig. [2]). 

The halo mass function depends sensitively on cosmo- 
logical parameters, including f2 m , and the equation- 
of-state of the dark energy. For our purposes, one of the 
most important cosmological parameters is the ampli- 
tude of the initial density fluctuations, characterized by 
as , the RMS variance of the linear density field, smoothed 
on scales of 8 Mpc. Uncertainty in this quantity trans- 
lates directly into uncertainty in the amplitude of the 
mass function. We utilize the latest value from WMAP, 
which provides a ~ 4% measurement of <7g j25| . For 
reference, a 5% error on erg shifts the contours in Fig- 
ure [T] by less than la in mass for a full-sky survey, and 
considerably less for smaller surveys. Since the value of 
cts is a major source of uncertainty in the use of the 
cluster mass function to constrain cosmology, there is 
great interest in improving its measurement. In addi- 
tion, the mass function also depends implicitly on the 
Hubble constant, h, which can be seen by expressing it 
in units of # of halos/(Mpc//i) 3 (observations naturally 
measure volume in these units). For simplicity we have 
explicitly put in the WMAP7 value (h = 0.710), but it is 
straightforward to re-express all of our results explicitly 
in terms of h (see Eqs.[T]and[2J and the text immediately 
beneath). 

The mass function predicts the number density of mas- 
sive dark matter halos in the Universe. For the purposes 
of this paper we are also interested in the scatter in this 
relation. At the high-mass end of the mass function, 
where the number density satisfies roughly one per vol- 
ume of interest, we assume that the distribution of ha- 
los is given by Poisson statistics. This is valid as the 
largest objects are spatially independent on these scales 
(> Gpc), and are dominated by shotnoise [26l[27|. 

We use the mass function presented in Tinker ct 
al. [lU, which gives the expected number density of dark 
matter halos, dn/dM, in units of Mpc/h, where h is the 
Hubble constant and volume is measured in comoving 
Mpc 3 . This mass function describes the abundance of 
sphcrical-overdensity dark matter halos, and is accurate 
to < 5% over the redshift range of interest (0 < z < 2), 
and for overdensity values (compared to the mean matter 
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FIG. 1: Contour plot of the most massive object in the Uni- 
verse. Three sets of contours are provided, for three differ- 
ent surveys: full sky, 178 deg 2 (corresponding to SPT), and 
lldeg 2 (corresponding to XMM2235) . The shaded contours 
represent the la and 2a (and for the lldeg 2 case, 3a) re- 
gions of the most massive halo in a ACDM Universe. The 
solid line contours are for the 2nd most massive halo, while 
the dashed line contours are for the 3rd most massive halo. 
The (blue) plus signs are Abell 2163 (double point) and Abell 
370, the three (green) diamonds are the three most massive 
clusters in the SPT 178 deg 2 survey, and the (red) square is 
XMM2235. Note that the mass values for Abell 2163 span the 
predicted region, while Abell 370 is slightly high. The SPT 
masses fit within their respective contours, while XMM2235 
is well outside its 2a contour. All masses are M200: spherical 
overdensity halos with A = 200 (measured with respect to 
Pmattcr). For data measured using different overdensities, we 
have converted to the M200 value which gives the equivalent 
probability. 

density at z) in the range 200 < A < 2300. This mass 
function has been calibrated for M200 £ 4 x 10 15 M Q , 
and therefore the extreme high-end of our calculations 
relies on extrapolation. In what follows we assume the 
WMAP7 cosmological parameters, namely, h = 0.710, 
O m = 0.264, n A = 0.734, and a s = 0.801 [||. 

The most massive object (Theory) — We are interested 
in determining the mass of the most massive object in 
our Universe. We calculate the expected distribution of 
masses at the high mass end, assuming Poisson statistics; 
the results are shown in Figure [1] The most massive ob- 
ject in the Universe is expected to be found at z = 0.22, 
with a mass M200 = 3.8 x 10 15 M Q . The marginalized 
la range in mass is 3.3 x 10 15 < A/200 < 4.4 x 10 15 , 
while in redshift it is 0.12 < z < 0.36. If the most 
massive object in the Universe falls outside the range 
2 x 1O 15 M < M 200 < 10 16 M o , we can conclude with 
high confidence that cither the initial conditions are non- 
Gaussian, or the growth of structure deviates from the 
predictions of general relativity. 

Figure [1] includes contours of the 2nd and 3rd most 
massive halos in the Universe. Going from the most mas- 
sive to the 2nd most massive results in a noticeable shift, 
demonstrating the power of just a few halos to constrain 
cosmology. As we go further down (e.g., from the 2nd 
to the 3rd most massive), the contours rapidly converge 
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due to the exponential steepening in expected number at 
lower mass. Note that the most massive halo occurs at 
low redshift. Furthermore, the contours are not centered 
on the most likely point; there is much larger scatter to 
high mass, with a sharp lower mass limit, due to the ex- 
ponential steepening. Note that these likelihoods are not 
independent, since if the most massive object has an un- 
usually low mass, it is assured that the subsequent few 
most massive objects will also be unusually low. We have 
performed Monte-Carlo studies which show that the cor- 
relations are weak, however, and the distribution of sepa- 
rations is well approximated by assuming the likelihoods 
are drawn independently. Figure [T] also shows the con- 
tours for the 1st and 2nd most massive objects from the 
recent SPT 178 deg 2 survey [HI, as weu as the contours 
for the archival XMM-Newton survey which discovered 
XMM2235. 

Figure [5] shows contours of the expected number of 
halos greater than a given mass, and found beyond a 
minimum redshift: (TV) (> M200, > z). The contours are 
roughly linear in the redshift range 0.2 < z < 2, and are 
well approximated (to better than 5%) by the family of 
lines: log 10 (M(A0 z)) = a{M) + b{M)z with 

a (TV) = 15.72 - 0.136A/"-0.014AA 2 -0.0012AA 3 (1) 
6(A0 = -0.5375 + 0.0058LV+ 0.0024A/' 2 + 0.00027A/' 3 , 

where N = log 10 (TV). For the redshift range z < 0.2, the 
results are well represented by the values at z = 0, which 
are given (to better than 2%) by: 

log 10 (M(A0) = 15 - 6 ~ 0.142AT - 0.014A/" 2 . (2) 

These expressions can be utilized to calculate the ex- 
pected number of objects above a given minimum mass 
and redshift in the mass range 1O 14 M < M200 < 
10 16 M Q and redshift range < z < 2, for any survey 
size. For a volume-limited sample, we are interested in 
(TV) (> M 2 oo, < z). These contours start at at z = 
(since there is no volume), and rapidly rise to their max- 
imum values, flattening by z ~ 0.2 at the values given 
by Eq. H Note that Eqs. [Q and M assume the WMAP7 
value of the Hubble constant, h = 0.710. To explicitly 
put in the h dependence, M200 and (TV) can be rescaled 
by (ft/0.71) and (0.71/h) 3 , respectively. 

The most massive object ( Observations) — The most mas- 
sive object in the Universe is likely to have already 
been detected by HQS AT (potentially even if it is be- 
hind the galactic plane (30[). Reliably measuring the 
masses of candidate ROSAT sources remains challenging, 
however, and therefore the specific identity and mass of 
the most massive object is unknown at present. Per- 
haps the most compelling candidate is Abell 2163 at 
z = 0.203, which has an X-ray mass measurement of 
Mgooc = 3.4 ± 0.8 x 10 15 M Q P, HI (where "500c" in- 
dicates A with respect to p cr a rather than p ma tter)- We 
expect 0.02 (0.002/0.2) clusters with at least this mass 
and redshift in the entire Universe, where the numbers in 
parentheses are the ler lower and upper bounds on (TV). 
An alternative, weak lensing measurement of the mass 
yields a lower value of M 50 oc = 2.0 ± 0.3 x 10 15 M® [H[, 
which has expectation 1.4 (0.5/4) (precisely agreeing 




0.0 0.5 1.0 1.5 2.0 



FIG. 2: Expected number of halos at redshift > z m - m with 
mass > M2oo,min, for a full sky survey. Each contour line 
represents a value of log 10 (TV). For a survey with fraction, /, 
of the full sky, the expected numbers of halos are diminished 
by the factor /. The dashed (red) line shows the result for 
(TV) = 0.01 using the fit from [20], based on an FOF halo 
finder with b = 0.2. It is virtually indistinguishable from the 
corresponding SO (A = 200) contour. The dotted (red) line 
represents the (TV) = 1 contour for a A = 200c mass function, 
with overdensity compared to p cr it, instead of the average 
matter density, p ma tter. Note that this agrees with the fiducial 
"0" line (A = 200) at high redshift, as the Universe becomes 
matter dominated. The data points are the same as in Fig. [T] 
Fitting forms for the curves in this figure are provided in the 
text. 



with predictions). Furthermore, [34| find an X-ray mass 
of M 5 ooc = 2.3 ± 0.07 x 10 15 M Q , which agrees well with 
the lensing value. Abell 370 is another compelling can- 
didate, with a weak lensing mass of M V i r = 2.93+q 32 x 
10 15 H~ 1 Mq at z = 0.375 [HI, [Hj], and an expectation of 
0.02 (0.005/0.05). These data points are shown in Fig- 
ures [T] and [H where we have converted the masses to the 
M200 values which give the equivalent probabilities. 

The figures also show the three most massive clusters 
from the SPT 178 deg 2 survey [28|, where we have added 
the statistical and systematic errors in quadrature. For 
the most massive cluster (M 200 = (8.3±1.7) x 10 14 M Q /h 
at z = 0.8), we would expect 0.14 (0.04/0.5) clusters in 
the given sky area with a mass and redshift at least as 
large. For the 2nd most massive cluster (A/200 = (8-2 ± 
1.9) x 10 14 M Q /Ti at z = 0.3), the expected number goes 
up to 2 (0.8/6), while for the 3rd most massive (M 2 oo = 
(6.56±1.54)xlO 14 M /fratz = 0.32) we expect 5 (2/14). 
These masses are fully consistent with theory. 

We also plot XMM2235, with a mass of M 20 o c = 
(7.3 ± 1.3) x 10 14 Mq at z = 1.4 Qj|- This cluster was 
found in an 11 deg 2 survey (/ = 0.0003). From Figured 
we would expect to find a few thousand objects with at 
least this mass in the entire Universe (z > 0), and only 10 
such objects at z > 1.4 on the entire sky. The expected 
number of clusters in an 11 deg 2 survey, with this mini- 
mum mass and redshift, islxl0~ 3 (3xlO~ 4 /4xlO~ 3 ). A 
conservative lower limit of M324 = 5 x 10 14 M Q is quoted 
in [l6|, which leads to an expectation of 6 x 10 -3 in the 
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survey area (see also [13, EH)- From Figure Q] we see that 
XMM2235 is a 3a outlier. Alternatively, the cluster's 
true mass would have to be reduced by 4a to achieve 
(N) = 1 (see Figure ^ . We note that these results are 
relatively insensitive to errors in the mass determination; 
15% errors do not qualitatively alter our conclusions. 

Current data argues for further exploration of the 
highest-mass end of the mass function, both at low and 
high redshift. It would be particularly difficult, the- 
oretically, to account for excessively massive clusters 
at z > 1, while having agreement at lower redshift 
(e.g., non-Gaussianity would not suffice). We expect to 
have dramatically improved complete high-redshift clus- 
ter surveys with which to test ACDM in the near future, 



including the full SPT survey (2000 deg 2 ), the Dark En- 
ergy Survey (5000 deg 2 ), Planck (all-sky), and eventually 
LSST (20, 000 deg 2 ). In particular, Planck is expected to 
provide a relatively complete, all-sky survey of all mas- 
sive clusters out to high redshift in the near future f37| . 
If the results from these cluster surveys disagree with the 
predictions outlined above, the ACDM paradigm for the 
growth of structure will need to be revisited. 
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Jerry Jungman and Michael Warren. 
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